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1. INTRODUCTION 

One of the fundamental questions that is open to both theory and experiment is 
how to observe the chiral phase transition. Indications that a phase transition 
from a chirally symmetric to a chirally broken phase have been obtained a 
long time ago: lattice gauge simulations of quantum chromodynamics (QCD), 
which are equilibrium calculations that can include only temperature but not 
finite density, show a phase transition at a finite critical temperature T c 
Unfortunately, due to the difficulties inherent in a confining theory such as 
QCD, it is impossible to observe this transition directly experimentally, let 
alone obtain experimental values for the usual quantities that are associated 
with critical phenomena, such as critical exponents. 

The experimental search for clues to a phase transition that has been con- 
ceived, is to proceed via heavy ion collisions, in which both high temperatures 
and densities are to be reached. Much effort has been expended in this direc- 
tion and will continue to be so during the early course of the next century, 
particularly at RHIC and CERN. However, from the theoretical side, there is 
at present no single conclusive observable that might indicate the occurrence 
of the chiral phase transition. In addition, the simple equilibrium picture of 
the static lattice simulations is inadequate for describing heavy ion collisions. 
Thus we may speculate as to whether a concrete signal of this phase transition 
can be obtained from a non-equilibrium formulation of the problem or trans- 
port theory, that is based on an underlying chiral Lagrangian. Here one may 
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surmise whether phenomena such as critical scattering, i.e. divergence of the 
scattering cross-sections at the phase transition, may leave an obvious signal, 
or simply the fact that the dynamically generated quark mass is a function of 
space and time may leave some visible remnant of a phase transition. 

Thus our goal is to construct a consistent transport theory that is based on 
the Nambu-Jona-Lasinio (NJL) model, which to order 1/N C also includes the 
mesonic degrees of freedom and the hadronization of quarks into these degrees 
of freedom. In this report, we discuss primarily the first term in the \/N c 
expansion in the collision, explicitly relating this to Feynman diagrams for 
quark-quark and quark-antiquark scattering (for details, see Ref.[||), leaving 
the mesonic sector to a future publication EJ. 



2. GENERAL TRANSPORT THEORY FOR FERMIONS. 

Non-equilibrium phenomena are completely described by the Schwinger - Kel- 
dysh formalism for Green functions Many conventions exist in the litera- 
ture. For our purposes, it is simplest to use the convention of Landau Q. In 
this, the designations + and — attributed to the closed time path that is shown 
in Fig.l, and the fermionic Green functions are defined as 

iS c (x,y) = (Tip(x)-tp(y)) = iS (x,y) 

iS a {x,y) = (fiP(x)iP(y))=iS++(x,y) 

iS > (x,y) = (i>(x)i>(y)) =iS + ~(x,y) 

iS<{x,y) = -{^{y)^{x))=iS- + {x,y). (1) 

In a standard fashion, one constructs equations of motion for the matrix of 
Green functions and one moves to relative and centre of mass variables, u = 
x — y and X = (x + y)/2. A Fourier transform with respect to the relative 
coordinate, or Wigner transform, 

S(X,p) = Jd\e^s(x+^X-^) (2) 

is then performed. Of particular interest is the equation of motion for S h = 
S' < . Regarding this, together with the equation of motion for the adjoint 
function S <<! , one arrives at the so-called transport and constraint equations 
by adding and subtracting these. One finds 

i{r,%^} + iAs- + (x,p)] = i- (3) 

and 

it) F)^ — ^ 
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respectively. Here 

h = Icon + I$ + 1%, (5) 

where the collision term is 

/con = ^- + {x, P )ks + -{x,p)-n + -{x,p)ks- + {x,p) 

_ Tgain r loss ( R \ 
— -'coll J coll > W 

and terms containing retarded (R) and advanced (A) components are contained 
in 

j| = -Z-+(X,p)AS R (X,p) ± S R (X,p)AX-+(X,p) (7) 

and 

I* = £ A (A»AS-+(X,p) T S-+(X,p)AZ A (X,p). (8) 
In these equations, A is given as 

Note that the equations (J3j) and fl) are general equations that are generic for 
any relativistic fermionic theory km. 



3. APPROXIMATION SCHEMES - APPLICATION TO THE NJL 
MODEL 

A double expansion in inverse powers of the number of colors l/N c and ft is 
performed. We explain this briefly in the context of the model. 
The simplest form of the NJL Lagrangian 

C = i>{ih ft-m )ip + G [(^0) 2 + (^7 5 ^) 2 j , (10) 

is considered. Here G is a coupling strength and mo is the current quark mass. 
The color degree of freedom is implicit. As this is a strong coupling theory, an 
expansion in G is inadmissable, and the usual technique involves an expansion 
in 1/N C . In this ordering, the lowest order Feynman graph is a Hartree term, 
which describes the mean field experienced by the quarks, while higher orders 
introduce collisions. We examine these two cases individually. In order to 
make connection with the classically known transport equations, we consider 
an expansion of the transport and constraint equations in powers of Ti. 

In the transport and constraint equations, Eqs.(|^) and ([|), the right hand 
sides I- and i+, respectively depend on h in two ways: (i) The derivative 
operator A, Eq. (|9]), contains Ti to all orders. Its expansion in % generates the 
various terms of the so-called gradient expansion, and (ii) the self-energies T, R 



R _ 



in. 



and S + may have overall factors in ft. While to order 0((1/N C ) ), E 
and is directly related to the mass, which we treat as a quantity that has a 
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direct classical interpretation, we find £ h oc hda/dft, and is only related via 
a factor Ti to a quantity with direct physical interpretation. Note that while 
we use the overall factors in h° or h 1 in or £ h in our fi classification, we 
do not expand the quantities m or da/dfl in powers of H. We will detail this 
strategy in the calculation of the mean field term and collision integrals. 

3.1. The Hartree Approximation 

The leading term in the 1/N C expansion corresponds to the Hartree approxi- 
mation. The Hartree self-energy can be evaluated, given a form for the Green 
function. Here we assume that 

S H +(X,p) = 2iri^i5(p -E p )Y,Us(p)u s (p)f q (X,p) 

S 

+S( PQ + E p ) Vs(-p)v s (-p)fdX, -P) (11) 

S 

with El — p 2 + m 2 (X) and f q ^ q the (spin independent) quark and antiquark 
distribution functions. This is the quasiparticle assumption. Similar expres- 
sions can be written down for the other Green functions of Eq.(l). Inserting 
the quasiparticle ansdzte into the Wigner transformed expression for the self- 
energy in the Hartree approximation, — i£^~ = 2i%ttiS^j~{x, x) leads to the 
non-equilibrium gap equation 

m{X) = m +4GN c m(X) J (2 ^ )3 E \ x) f 1 ~ M X >P) ~ fi(X,P% ( 12 ) 

Some manipulations are now necessary for the evaluation of (3) and (4). The 
evaluation of I± is made according to the % expansion. Then, on performing 
a spinor trace and integrating over positive or negative energies, leads one to 
the Vlasov equation 

p^d^j^X.p) + m(X)d fM m(X)d^f qtq (X,p) = 0, (13) 

where p° = E p (X), and the constraint equation 

(p 2 -m 2 (X))f q ^ = Q. (14) 

to leading order in Ti. Note that the constraint equation validates the quasi- 
particle ansatz, and that the equations Eq. (|l3|) and jl^) can be solved self- 
consistently. 

3.2. The Collision Term 

There are several graphs which contribute to the next order in the 1/N C ex- 
pansion. One of these is displayed in Fig. 1(a). The accompanying graph of 
Fig. 1(b) is of higher order in l/N c , and is not required on this basis, but, as 
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Fig. 1. — Contributions to the self-energy that are of higher order in 1/N C than the 
Hartree term 



it turns out, is essential for a complete and correct identification of Feynman 
amplitudes in casting the collision term into a Boltzmann-like form. The ini- 
tial task is the evaluation of 7 co n of Eq.(6), with scalar or pseudoscalar vertices. 
There are three possible combinations, T, a , and E C oii, which contain unity 
and «75 throughout or unity or 175 in a mixed fashion. For example, the scalar 
self-energy is determined to be 

^+-fv \ A^2f,2 f d4 Pi d4 P2 d 4 p 3 4 

x [S + -(X, Pl )tr (s- + (x, P2 )s + -(x,p 3 )) 

-S + -(X, Pl )S-+(X, P2 )S+-(X, P3 )], (15) 

while the other self-energies are similar 

As before, performing a spinor trace and integrating over a positive or nega- 
tive energy range is required, so that the loss term, for example, is correspond- 
ingly given as 



T «,l0SS 



2ni 



'coll 



2E p 



'(X, Po = E p ,p)^2u s ( P )u s (p)f q (X,p) 



(16) 



which, when evaluated using the quasiparticle approximations for the Green 
functions, leads to eight terms. After some algebra, one can show that only 
terms that represent elastic scattering survive energy integration. Other terms 
that contain vacuum fluctuations, as well as particle production and annihila- 
tion are present at this level, see Fig. 2, but are suppressed by the quasiparticle 
approximation. An off-shell ansatz is required for them not to vanish. It is also 
important to note that both diagrams of Fig. 1 are essential for the identifica- 
tion of the Feynman amplitudes for elastic scattering in the s, t and u channels, 
as appropriate. While graph (a) leads to the square of the individual ampli- 
tudes, graph (b) is necessary for the mixed terms that allow one to construct 



6 



S.P. Klevansky et al. 



Fig. 2. — Processes of the form q — ► (qq)q and (qg) — > q, as well as vacuum fluctu- 
ations that are suppressed in the quasiparticle approximation. These are heuristic 
graphs and are not Feynman graphs 



a matrix element squared that involves several channels. The net result is 

p»dMX,p) +m(X)d^m(X)dZf q (X,p) = 

x {^|„- OT (p2 - mfM)fMfM)JM - f(pi)fMfMf q (p)) 
+^l M -- M -(p2 - imMJMfMW) - f q ( P i)fMfMf q (p))}, 

(17) 

where the scattering cross-sections contain s, t and u channel graphs as appro- 
priate, in the Born approximation and which is the relativistic generalization 
of the result of Kadanoff and Baym. To lowest order in K, arguments can be 
given to suggest that the constraint equation remains unaltered. 

4. DISCUSSION 



The most important facets of this calculation arc (i) in identifying precise 
Feynman amplitudes that contribute to the scattering processes in the collision 
term, and which represents the relativistic generalization of the Kadanoff-Baym 
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work [0 , and (ii) the recognition of additional processes that would be present 
should a quasiparticle assumption be relaxed, such as vacuum fluctuations and 
particle production and annihilation. The double expansion in Ti and l/N c 
lays the foundation for including additional higher order terms in the 1/N C 
expansion, from which mesons are constructed, and which enables one to build 
a dynamical model for both mesons and quarks. For this, we refer the interested 
reader to §, §. 

We have started numerical simulations of this chiral transport theory. A first 
study that solves the Vlasov equation by means of a finite difference technique 
and which inserts collision terms in the relaxation time approximation for three 
flavors is presented in Ref. Q. In this, one finds that the collisions determine 
the high energy parts of the particle distribution, which is exponential. The 
mean field, on the other hand, produces an enhancement over the exponential 
fits at low energies. The behavior of the quark mass as a function of time is 
such that a dip occurs in the initial state, and this gradually flattens as the 
quarks flow out. 
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